home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2002 November / SGI IRIX Base Documentation 2002 November.iso / usr / share / catman / p_man / cat3 / SCSL / dlagv2.z / dlagv2
Encoding:
Text File  |  2002-10-03  |  4.6 KB  |  133 lines
  1.  
  2.  
  3.  
  4. DDDDLLLLAAAAGGGGVVVV2222((((3333SSSS))))                                                          DDDDLLLLAAAAGGGGVVVV2222((((3333SSSS))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DLAGV2 - compute the Generalized Schur factorization of a real 2-by-2
  10.      matrix pencil (A,B) where B is upper triangular
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE DLAGV2( A, LDA, B, LDB, ALPHAR, ALPHAI, BETA, CSL, SNL, CSR,
  14.                         SNR )
  15.  
  16.          INTEGER        LDA, LDB
  17.  
  18.          DOUBLE         PRECISION CSL, CSR, SNL, SNR
  19.  
  20.          DOUBLE         PRECISION A( LDA, * ), ALPHAI( 2 ), ALPHAR( 2 ), B(
  21.                         LDB, * ), BETA( 2 )
  22.  
  23. IIIIMMMMPPPPLLLLEEEEMMMMEEEENNNNTTTTAAAATTTTIIIIOOOONNNN
  24.      These routines are part of the SCSL Scientific Library and can be loaded
  25.      using either the -lscs or the -lscs_mp option.  The -lscs_mp option
  26.      directs the linker to use the multi-processor version of the library.
  27.  
  28.      When linking to SCSL with -lscs or -lscs_mp, the default integer size is
  29.      4 bytes (32 bits). Another version of SCSL is available in which integers
  30.      are 8 bytes (64 bits).  This version allows the user access to larger
  31.      memory sizes and helps when porting legacy Cray codes.  It can be loaded
  32.      by using the -lscs_i8 option or the -lscs_i8_mp option. A program may use
  33.      only one of the two versions; 4-byte integer and 8-byte integer library
  34.      calls cannot be mixed.
  35.  
  36. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  37.      DLAGV2 computes the Generalized Schur factorization of a real 2-by-2
  38.      matrix pencil (A,B) where B is upper triangular. This routine computes
  39.      orthogonal (rotation) matrices given by CSL, SNL and CSR, SNR such that
  40.  
  41.      1) if the pencil (A,B) has two real eigenvalues (include 0/0 or 1/0
  42.         types), then
  43.  
  44.         [ a11 a12 ] := [  CSL  SNL ] [ a11 a12 ] [  CSR -SNR ]
  45.         [  0  a22 ]    [ -SNL  CSL ] [ a21 a22 ] [  SNR  CSR ]
  46.  
  47.         [ b11 b12 ] := [  CSL  SNL ] [ b11 b12 ] [  CSR -SNR ]
  48.         [  0  b22 ]    [ -SNL  CSL ] [  0  b22 ] [  SNR  CSR ],
  49.  
  50.      2) if the pencil (A,B) has a pair of complex conjugate eigenvalues,
  51.         then
  52.  
  53.         [ a11 a12 ] := [  CSL  SNL ] [ a11 a12 ] [  CSR -SNR ]
  54.         [ a21 a22 ]    [ -SNL  CSL ] [ a21 a22 ] [  SNR  CSR ]
  55.  
  56.         [ b11  0  ] := [  CSL  SNL ] [ b11 b12 ] [  CSR -SNR ]
  57.         [  0  b22 ]    [ -SNL  CSL ] [  0  b22 ] [  SNR  CSR ]
  58.  
  59.         where b11 >= b22 > 0.
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDLLLLAAAAGGGGVVVV2222((((3333SSSS))))                                                          DDDDLLLLAAAAGGGGVVVV2222((((3333SSSS))))
  71.  
  72.  
  73.  
  74. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  75.      A       (input/output) DOUBLE PRECISION array, dimension (LDA, 2)
  76.              On entry, the 2 x 2 matrix A.  On exit, A is overwritten by the
  77.              ``A-part'' of the generalized Schur form.
  78.  
  79.      LDA     (input) INTEGER
  80.              THe leading dimension of the array A.  LDA >= 2.
  81.  
  82.      B       (input/output) DOUBLE PRECISION array, dimension (LDB, 2)
  83.              On entry, the upper triangular 2 x 2 matrix B.  On exit, B is
  84.              overwritten by the ``B-part'' of the generalized Schur form.
  85.  
  86.      LDB     (input) INTEGER
  87.              THe leading dimension of the array B.  LDB >= 2.
  88.  
  89.      ALPHAR  (output) DOUBLE PRECISION array, dimension (2)
  90.              ALPHAI  (output) DOUBLE PRECISION array, dimension (2) BETA
  91.              (output) DOUBLE PRECISION array, dimension (2)
  92.              (ALPHAR(k)+i*ALPHAI(k))/BETA(k) are the eigenvalues of the pencil
  93.              (A,B), k=1,2, i = sqrt(-1).  Note that BETA(k) may be zero.
  94.  
  95.      CSL     (output) DOUBLE PRECISION
  96.              The cosine of the left rotation matrix.
  97.  
  98.      SNL     (output) DOUBLE PRECISION
  99.              The sine of the left rotation matrix.
  100.  
  101.      CSR     (output) DOUBLE PRECISION
  102.              The cosine of the right rotation matrix.
  103.  
  104.      SNR     (output) DOUBLE PRECISION
  105.              The sine of the right rotation matrix.
  106.  
  107. FFFFUUUURRRRTTTTHHHHEEEERRRR DDDDEEEETTTTAAAAIIIILLLLSSSS
  108.      Based on contributions by
  109.         Mark Fahey, Department of Mathematics, Univ. of Kentucky, USA
  110.  
  111.  
  112. SSSSEEEEEEEE AAAALLLLSSSSOOOO
  113.      INTRO_LAPACK(3S), INTRO_SCSL(3S)
  114.  
  115.      This man page is available only online.
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.